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ABSTRACT 

We report that the Sloan Digital Sky Survey quasar SDSS J165043.44+425149.3 is 
gravitationally lensed into two images, based on observations obtained with the WIYN 
3.5 m telescope at the Kitt Peak National Observatory. The lensed quasar, at a redshift 
of z = 1.54, appears as two images separated by 1'.'2 with 5-band magnitudes of 17.8 
and 20.0. The lensing galaxy is clearly detected in /-band images obtained in / /3 
seeing, after point spread function subtraction of the two quasar images. A strong 
metal-line absorption system is also identified in the unresolved SDSS spectrum of the 
double quasar, suggesting a plausible lens redshift of z = 0.58. The UBRI flux ratios 
of the pair vary significantly from 8.5:1 in the blue to 5.4:1 in the red, a difference of 
0.5 magnitudes, and may indicate the presence of microlensing in one or both quasar 
components. The predicted differential time delay between quasar images is on the 
order of one month, assuming the intervening absorption system is due to the lensing 
galaxy. 

Subject headings: gravitational lensing: individual (SDSS J165043. 44+425149. 3) 



1. Introduction 

Gravitationally lensed quasars are useful for a wide range of applications, such as constraining 
the cosmological parameters H Q (Koopmans et al. 2000; Kochanek et al. 2002) and A (Falco et al. 
1998; Chae et al. 2002), mapping the structure of luminous and dark matter in intermediate redshift 
galaxies (Koopmans & Treu 2003), and studying the evolution of the lens galaxies themselves 
(Kochanek et al. 2000). Owing to a number of large-scale surveys for lensed quasars at both optical 
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and radio wavelengths, there are now well over 70 confirmed systems. In this paper, we report that 
the quasar SDSS J1650+4251 is gravitationally lensed into two images, based on observations 
obtained with the WIYN 3.5 m telescope. We present the lens discovery data and subsequent 
detection of the lensing galaxy in §2, assess the SDSS spectrum of the system in §3, and discuss 
the plausible redshift and potential model of the lensing galaxy in §4. Our findings and possible 
directions for future work are summarized in §5. 

2. Observations and Analysis 

To the best of our knowledge, 5 SDSS J1650+4251 (16 h 50 m 43!5, +42° 51' 4570; J2000.0) was 
first identified as a z = 1.54 redshift quasar, with a blue magnitude of g = 17.6, from the first 
data release of the Sloan Digital Sky Survey (SDSS DR1; http://www.sdss.org/drl). During 
the nights of 24 and 25 April 2003, approximately 200 DR1 quasars were reimaged under excellent 
seeing conditions with the WIYN telescope at the Kitt Peak National Observatory, with the goal of 
identifying multiple image morphology indicative of gravitational lensing. Our snapshots consisted 
of 30 second i?-band exposures using the Mini-Mosaic direct imager, which has a field of view of 
9f 6 square and a gain and readnoise of 1.5 e~ ADU -1 and 5.5 e~, respectively. We operated in 
2x2 binning mode to help minimize the readout time, which provided a scale of / /282 pixel -1 . 

For each target, we assigned an a priori lensing probability based on the quasar's redshift 
and g-band magnitude. These crude probabilities, which are computed using an optical depth 
calculation and a magnification bias factor (e.g., Kochanek 1996), were used to select targets for 
reimaging during the observing run. The a priori lensing probability for SDSS J1650+4251 was 
0.7%, placing it in the top 2% of the ~8000 DR1 quasars accessible during the WIYN run. The 
discovery snapshot of the target immediately revealed it to be double with an image separation 
of l'/2 and an i?-band flux ratio of ~ 5 : 1. Binned followup images were obtained in Harris 
U BRI filters, as well as several unbinned images in /-band alone, and revealed a similar separation 
and flux ratio in the other bandpasses. The seeing during the discovery and followup images was 
exceptional, ranging from 076 FWHM in U to 073 FWHM in /. A summary of the WIYN data 
for SDSS J1650+4251 is provided in Table 1. In addition, UBRI images of the Landolt (1992) 
standard field PG 1633 were also obtained immediately following the quasar observations. A 2'.5 
square subraster of the target and surrounding field is shown in Figure 1. 

The images were bias-subtracted and flat-field corrected using twilight sky flats taken during 
the run. To model the light distribution of the double image, we have simultaneously fit two 
empirical point spread functions (PSFs) to components A (brighter) and B (fainter), for each 
image listed in Table 1, using a Powell minimization routine (Press et al. 1992). For the BRI-b&nd 



5 Queries to the NRAO VLA Sky Survey radio catalog, the ROSAT X-ray database, and the NASA/IPAC Extra- 
galactic Database returned no results within 2' of the quasar's coordinates. 
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images, star #2 in Figure 1 provided the PSF. This star dropped significantly in brightness toward 
the blue, so star #1 served as the [/-band PSF. After subtracting the best-fit models, no significant 
structure was detected in the residual image for the UBR data. The average flux ratio for the 
two components was 5.4:1 in R, and increased toward bluer wavelengths to 7.3:1 in B and 8.5:1 
in U. Relative astrometry and photometry for the two-component models are reported in Table 2. 
Magnitudes in Table 2 have been calibrated using zeropoints (but no color terms) obtained from 
the PG 1633 standards. 

The lack of significant residuals in the U BR data indicates that the system is well-modeled by 
two point sources at these wavelengths. However, significant structure was present at the position 
of the double image after subtracting the two-component model from the unbinned /-band data, as 
seen in panel (b) of Figure 2. The peak of the residual flux is ~15o" (where a is the rms noise per 
pixel), and is peaked roughly collinear and between the brighter two components. This position is 
what would be expected for a simple singular isothermal sphere (SIS) lens model of the system, 
assuming the two brighter and point-like images were lensed components of a single background 
quasar and the faint residual flux was from a foreground lensing galaxy. 

To model the suspected galaxy flux (hereafter, component G), a third component was added 
to the /-band models and the simultaneous fits were repeated. The galaxy was modeled using a 
circularly-symmetric "pseudogaussian" light profile (Schechter & Moore 1993), which falls off slower 
than a true gaussian at moderate distances (several pixels) from the profile core. The relative fluxes 
and positions of the three components were free to vary. Panel (c) of Figure 2 shows the stacked 
residuals from the three-component /-band fits, and is free from significant structure above 5a at 
the position of component G. There are still significant residuals (±8cr) at the cores of component 
A and B, but this is most likely an artifact of fractional pixel shifts of the undersampled PSF. 
Panel (d) shows the same fit as in panel (c), but with the galaxy left unsubtracted. Results for the 
three-component fit are also reported in Table 2. 

The /-band data place the galaxy closer to component B (6*bg=0 / . / 36) than to A (#ag=0 / '87), 
and offset by C//15 (1 pixel) to the south-west with respect to the quasar AB line. The average /- 
band flux ratio, 5.75:1, is slightly larger than that found at R. To estimate the galaxy magnitude in 
/, we have performed aperture photometry on the unsubtracted galaxy flux for the three-component 
model, that is, on panel (d) of Figure 3. An aperture radius of (X'71 (5 pixels) was used, and yields 
a calibrated magnitude of / = 20.5. Since no aperture correction was applied, this is likely a lower 
limit to the true galaxy magnitude. 

3. Spectral Observations 

Unresolved spectra of SDSS J1650+4251 was obtained with the 2.5 m telescope at Apache 
Point, New Mexico, on 19 June 2001 as part of the standard SDSS spectroscopic followup of color- 
selected quasar candidates (Richards et al. 2002). The spectrum has a wavelength coverage of 
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approximately 3800 A to 9200 A, a dispersion that ranged from 1 to 2 A pixel" 1 from blue to red, 
and a total on-source exposure time of 2700 seconds. A 5 A binned sample of the Sloan spectrum 
is reproduced in the top panel of Figure 3. 

The 50% seeing disc was ~ 2" during the exposure, ensuring that light from both components 
were well mixed. The unresolved spectrum is consistent with at least one object being a z = 
1.541 quasar, based on gaussian fits to the CIII] and Mgll broad emission lines. If the remaining 
component were a foreground star, then one might expect zero-redshift stellar features in the 
unresolved spectrum. The expected locations of several stellar features are marked by the dashed 
lines in Figure 3, but none are readily identified with features in the spectrum. (The absorption 
feature at the expected location of the Na doublet is actually a singlet). The lack of significant 
stellar absorption lines, together with the absence of other quasar emission lines at a redshift 
different than that quoted above, argues for two similarly redshifted quasars. 

There is also evidence for two z > absorption systems in the spectrum, based on four 
doublet features at approximate observed wavelengths of AA3945, 4090, 4415, and 7125 A. The 
bluest and reddest absorption doublets lie on top of the CIV and Mgll quasar emission lines, and 
are consistent with CIV AA1548,1550 (bottom left panel of Figure 2) and Mgll AA2796,2803 at 
redshifts of z = 1.545 and z = 1.546, respectively. These features are clearly associated with 
the quasar. The remaining two doublets are consistent with an intervening absorption system 
at z = 0.577 arising from Fell AA2587,2600 and Mgll AA 2796,2803 (bottom right of Figure 2). 
Assuming this absorber redshift, a weak Mgl A 2853 absorption singlet was also identified at the 
same redshift. 



4. Discussion 

4.1. Plausible Redshift of the Lensing Galaxy 

Although the lensing galaxy currently lacks a spectroscopic redshift, a plausible redshift of 
z = 0.577 is suggested by the intervening Mgll and Fell absorption system present in the quasar 
spectrum. One can statistically estimate the redshift of the lensing galaxy following the approach 
of Kochanek (1992). Using the image separation of 1'.'2 and quasar redshift of z = 1.541, we find a 
median lens redshift of z = 0.59 with a la interval of 0.36 < z < 0.82 (assuming matter and energy 
densities of 0.3 and 0.7 of the closure density, respectively). 

An additional estimate of the galaxy redshift can be obtained from the observed /-band mag- 
nitude. For a given redshift and assuming an isothermal sphere model for the lensing potential, the 
galaxy's velocity dispersion a is related to the observed image separation A9 by 

* = (!) 

J-^os ° 

(Narayan & Bartelmann 1999), where Di s and D os are angular diameter distances from the lens to 
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the source, and from the observer to the source, respectively. The galaxy's 5-band luminosity can 
be estimated from its velocity dispersion using the Faber- Jackson relationship, L/L* = (<t/ct*) 7 , 
where L* corresponds to a i?-band magnitude of M* = —19.7 + 5 log h. and we adopt a* = 225 km 
s _1 and 7 = 4.0 appropriate for elliptical galaxies. Finally, the galaxy's apparent magnitude m is 
given by 

™ A B(A bs) - M AB (A rcst ) = 5 log + 7.51og(l + z), (2) 

where D a i is the angular diameter distance from the observer to the lens. For calculating MAB(Kest), 
a spectral energy distribution (SED) for an early- type galaxy was obtained from S. Lilly (1997, 
priv. comm.), which consisted of interpolated and extrapolated values of the SEDs presented by 
Coleman, Wu, & Weedman (1980). The SED is then normalized to the Faber- Jackson luminosity 
at 4400(1 + z) A, which yields the predicted AB magnitudes. To transform back to standard 
magnitudes, we use / = Iab — 0.456 (Fukugita et al. 1995). 

Figure 4 shows the predicted galaxy /-band magnitude as a function of redshift (curved solid 
line), along with the la error in the predicted magnitude from the observed dispersion in the Faber- 
Jackson relationship (curved dashed lines; Dressier et al. 1987). The galaxy appears systematically 
fainter as it is placed at higher redshifts, up to z ~ 1.2, since the velocity dispersion required to 
reproduce the observed image separation becomes systematically smaller. Beyond z ~ 1.2, the 
required velocity dispersion (and therefore apparent magnitude) increases rapidly. At the observed 
galaxy magnitude oil = 20.5 (horizontal line), the formally best galaxy redshift is at either z = 0.49 
or z = 1.39. The statistical constraint on the lensing galaxy redshift (shaded region) supports a 
redshift under unity, as does the observed metal-line absorption system (vertical tick). Taken 
together, the data favors a lensing galaxy redshift in the range of 0.4 < z < 0.8, consistent with 
the absorption system at z = 0.577. 



4.2. Lens Models 

We have fit a lens model to the system, which is straightforward to model given the simple 
image geometry. The misalignment between components A, B, and G implies that some asymmetry 
must be present in the lensing potential. Our basic model is an SIS potential for the lensing galaxy, 
plus an external shear to break the circular degeneracy. It has a two-dimensional effective lensing 
potential (f> given by 

cj) = br- |r 2 cos2(#-# 7 ). (3) 

The major axis of the isopotential lines point along # 7 (measured E of N), which is also the direction 
(modulo 180°) of the perturbing mass responsible for the shear. Note that the sign convention for 
the shear term is opposite to that adopted by Wisotzki et al. (2002), whose SIS+shear model 
predicts a perturbing mass at ±90° from T In the limit of vanishing shear (7 — ► 0), the strength 
of the potential, b, corresponds to the Einstein radius of the system, which is also half of the observed 
image separation. There are five model parameters (including two for the unknown quasar source 
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position), and five constraints (the relative positions of components A and B with respect to the 
lensing galaxy and the A/B flux ratio). Thus there are no degrees of freedom and our best fit is 
characterized by % 2 = 0. 

Parameter values were obtained by minimizing residuals in the source plane. For a given 
parameter set, each quasar image is traced back from the image plane to the source plane using the 
lens equation, and the optimal model parameters are then found by minimizing residuals between 
the model source position and the backward projection of the image positions. 

Using the I-band positions and flux ratio as constraints, our best-fit model parameters are 
b = 0"665, 7 = 0.152, and 6> 7 = 136.8°, and a lensed source position of (Aa, A<5) = (-tf'136, tf'218) 
with respect to the lensing galaxy. The model predicts magnification factors of 4.94 and -0.86 for 
components A and B, where the negative sign denotes a parity flip, and a differential time delay 
between quasar images of 33.8 days for a Hubble constant of H Q = 75 km s _1 and assumed galaxy 
redshift of z = 0.577. At this redshift, the corresponding velocity dispersion for the lensing galaxy 
would be 185 km s _1 for an elliptical galaxy. The Faber- Jackson relationship predicts a l?-band 
galaxy luminosity of roughly 0.5 L*, which is not unreasonable. 

The shear strength is fairly large for a double system. If the source of shear was from a 
neighboring galaxy with the same velocity dispersion and redshift as the main lensing galaxy, then 
it would be located 2 / .'2 away. Indeed, there is a possible companion galaxy 3"0 south from the 
lensing galaxy (object G2 in Figure 2), however the position angle of the nearby galaxy (171.5° 
E of N) differs by 35° degrees from the best-fit shear direction. While this may indicate that the 
companion is actually a chance projection, it is hard to dismiss the notion that the nearby galaxy 
and large shear are associated. 

In other lensed quasars with large shears, a nearby group or cluster of galaxies at moderate 
distances (10" to 40") sometimes acts as the source of the shear (Kundic et al. 1997; Schechter et 
al. 1997; Kneib et al. 2000). For SDSS J1650+4251, there is one faint object 8" from the lensing 
galaxy at a PA of 110°, but the light profile is more consistent with a point source than a galaxy. 
Otherwise, there are no obvious galaxy concentrations near the shear position angle. 

5. Conclusions 

We have reported that the z = 1.54 quasar SDSS J1650+4251 is gravitationally lensed into a 
l"2 double. The lensing galaxy is clearly detected in /-band after PSF subtraction of the two quasar 
images, and is located closer to the fainter quasar image, as expected for a simple SIS model of the 
lensing potential. A strong metal-line absorption system is detected in the unresolved spectrum of 
the quasar pair, suggesting a plausible lens redshift of z = 0.577. Taking this as the lensing galaxy 
redshift, an SIS plus external shear model predicts a differential time delay of 33.8 days (H Q = 75 
km s _1 ). Resolved spectra of the lensing galaxy will be needed to obtain a confident galaxy redshift 
and time delay prediction. 
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The flux ratios for the system show a significant change (0.5 magnitudes) with wavelength, 
ranging from 5.4:1 in R to 7.3:1 in B to 8.5:1 in U. Since gravitational lensing is achromatic, one 
would ideally expect the flux ratio to be constant, but in practice the individual magnifications are 
affected by extinction and microlensing from dust and stars in the lensing galaxy. Since the flux 
ratio for SDSS J1650+4251 increases towards the blue, this does indeed suggest that perhaps one 
of these scenarios is correct. For at least one doubly lensed quasar, SBS 0909+532 (Kochanek et al. 
1997), the flux ratio is known to vary by up to a full magnitude across optical wavelengths, so the 
magnitude of variation observed for SDSS J1650+4251 is not unprecedented. Optical monitoring 
of the double quasar, particularly at blue wavelengths, will help determine if microlensing is indeed 
significant and if the lensed quasar is sufficiently variable for time delay studies. 

The authors would like to extend their appreciation to George Will (Kitt Peak night assistant) 
and Charles Corson (WIYN site engineer) for their help with the WIYN snapshot survey. We 
would also like to thank Charles Bailyn for comments on the manuscript. 
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Fig. 1. — 2'5 square finding chart for SDSS J1650+4251 and surroundings. The lens and two PSF 
stars are labelled. 

Fig. 2.— Subrasters of the SDSS J1650+4251 taken in 073 seeing with WIYN (scale of 07141 
pixel -1 ). The orientation and angular size of each panel are identical. Panel (a): stacked /-band 
image of the lens. Panel (b): residuals of panel (a) after fitting and subtracting two empirical PSFs 
to components A and B. Panel (c): residuals of panel (a) after fitting and subtracting two empirical 
PSFs plus circularly-symmetric pseudogaussian profile to model the lensing galaxy (component G). 
Panel (d): The same as panel (c), but with the galaxy model unsubtracted. The contrast level 
for panels (b), (c), and (d) range from -3a to 10<7, where a is the rms noise (Poisson + readout) 
associated with each pixel. Note the possible companion galaxy (G2) 3" south of the main lensing 
galaxy 

Fig. 3. — Top panel: archive SDSS spectrum of SDSS J1650+4251. The presence of two metal-line 
absorption systems, one at z = 0.577 and the other close to the quasar redshift at z = 1.545, 
are indicated by the tick marks. Dashed lines mark the expected location of zero-redshift stellar 
features. Bottom-left panel: Closeup of the CIV absorption doublet at the quasar redshift. Bottom- 
right panel: Closeup of the intervening Mgll absorption doublet at z = 0.577. 

Fig. 4. — Predicted /-band magnitude of the lensing galaxy as a function of the lens galaxy redshift 
(curved solid line), along with the la error range (curved dashed lines) expected from the scatter 
in the Faber-Jackson relationship. The horizontal line marks the observed /-band magnitude of the 
lensing galaxy (/ = 20.5), and the vertical tick marks the location of the intervening absorption 
system [z = 0.577). The shaded region (0.36 < z < 0.82) is the la statistical redshift range for the 
lensing galaxy 
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Table 1. Log of Observations for SDSS 1650+4251 



Frame 


Filter 


Exp. (s) 


FWHM (") 


Bin Mode 


302 


R 


30 


0741 


2 


304 


U 


60 


0"58 


2 


305 


u 


90 


tf'60 


2 


306 


I 


60 


0"38 


2 


307 


B 


60 


0"45 


2 


308 


I 


60 


0"32 


1 


309 


I 


150 


0731 


1 


310 


I 


150 


tf'32 


1 


311 


I 


150 


tf'33 


1 


313 


R 


60 


tf'40 


2 


314 


B 


60 


tf'47 


2 



Note. - Binning modes 1 & 2 yield a pixel scale of 
07141 pixel -1 k tf'282 pixel -1 , respectively. 
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Table 2. Astrometry and Photometry for SDSS 1620+4251 



B G 



Filter m A m B m G AR.A. (") ADec. (") AR.A. (") ADec. (") 

U 17.10 19.42 

B 17.80 19.96 

R 17.44 19.26 

I 17.15 19.05 20.5 



Note. - Harris magnitudes for components A, B, and G are denoted by itia, rnB, and ti%g, 
respectively. Relative astrometry for components B and G are with respect to component A. Error 
bars for the /-band positions are from the rms dispersion among the four unbinned /-band frames. 



0.212 -1.138 

0.212 -1.145 

0.223 -1.147 

0.223 ± 0.002 -1.163 ± 0.001 0.017 ± 0.032 -0.872 ± 0.026 
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